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Abstract 

Hierarchical clustering of dark matter halos is thought to describe well the large scale 
structure of the universe. The baryonic component of the halos is shock heated to the 
virial temperature while a small fraction of the energy flux through the shocks may 
be energized through the first order Fermi process to relativistic energy per particle. 
It has been proposed that the electrons accelerated in this way may upscatter the 
photons of the universal microwave background to gamma ray energies and indeed 
generate a diffuse background of gamma rays that compares well to the observations. 
In this paper we calculate the spectra of the particles accelerated at the merger 
shocks and re-evaluate the contribution of structure formation to the extragalactic 
diffuse gamma ray background (EDGRB), concluding that this contribution adds 
up to at most 10% of the observed EDGRB. 
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1 Introduction 


EGRET observations [1] showed that the universe is permeated by a back¬ 
ground of gamma radiation that seems to exceed the flux of gamma rays 
expected from cosmic ray interactions in our own Galaxy, as calculated using 
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theoretical models of the origin and propagation of cosmic rays. This excess 
has long been considered of extragalactic origin, and innumerable attempts to 
explain it in terms of different kind of sources have been made. Whether this 
radiation is the result of discrete unresolved extragalactic sources or rather a 
truly diffuse background is still unknown, and is matter of investigation for 
future gamma ray telescopes such as GLAST. 

It is somehow disturbing that the extragalactic origin for this background has 
been inferred from a combination of measurements and theoretical modelling 
of the diffuse galactic gamma radiation. In fact some authors [2] have proposed 
that observations can be also explained as a result of a population of galactic 
relativistic electrons upscattering the microwave and starlight radiations to 
gamma ray energies through inverse Compton scattering (ICS). These elec¬ 
trons would not be correctly accounted for in standard models of cosmic ray 
propagation in the Galaxy. 

If this radiation is in fact mainly extragalactic, its source or sources need to 
be found. While it is believed that blazars may contribute a large fraction of 
the extragalactic diffuse gamma ray background [4,5,3,6] (EDGRB), it is not 
yet clear whether they can saturate it (see, for example [7,8]). Around 1 GeV 
a non negligible contribution to the diffuse background might also come from 
normal galaxies [9]. 

In the last few years, clusters of galaxies have been proposed as sources of 
high energy gamma rays and in fact even as sources of the EDGRB. The first 
paper to propose this possibility is Ref. [10]. Some problems were identified in 
these calculations and discussed in Refs. [11,12]. A detailed calculation of the 
EDGRB due to clusters of galaxies was carried out in [13], where the authors 
concluded that not more than a few percent of the observed background of 
gamma rays could be accounted for in terms of hadronic interactions in clusters 
of galaxies. 

More recently, the authors in Ref. [14,15] have reproposed a connection be¬ 
tween the EDGRB and clusters of galaxies. More correctly the connection 
should exist between the EDGRB and the process of hierarchical large scale 
structure formation. The claim is that the whole EDGRB can be explained in 
terms of IGS of electrons accelerated at shocks during structure formation up 
to ultrarelativistic energies. Shocks form naturally during the merger of two 
halos that generate a new bigger structure. 

In a recent paper [16] we have studied in detail the process of acceleration and 
reacceleration of particles at shocks during mergers of clusters of galaxies, and 
we have proposed a semi-analytical method to evaluate in a self-consistent 
manner the Mach numbers of the shocks developed at each merger event. 
The Mach numbers are related in a unique way to the strength of the shocks 
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and therefore to their ability to accelerate particles. We adopt here the same 
method introduced there and apply it to the calculation of the contribution 
of structure formation to the EDGRB. 

Our conclusion is that at most 10% of the observed EDGRB can be explained 
by invoking the process described above. We discuss in detail the reasons for 
the difference between our results and those in [14,15]. 

The paper is planned as follows: in §2 we describe the basics of hierarchical 
clustering and the formation of shock surfaces during structure formation. In 
§3 we summarize the physics of shock acceleration and specialize the discussion 
to the case of merger shocks. In §4 we illustrate our calculations of the diffuse 
gamma ray background from structure formation. We conclude in §5. 


2 Structure formation and related shocks 


The standard theory of structure formation predicts that larger structures re¬ 
sult from the mergers of smaller structures, which on average are formed at 
earlier times. Press and Shechter [17] (hereafter PS) were the hrst to propose 
an efficient analytical description of the hierarchical clustering. It represents 
an extremely powerful tool that allows one to reconstruct realizations of the 
merger history of a cluster with hxed mass at the present time. There are now 
different flavors of these analytical methods with different levels of sophistica¬ 
tion [18,19]. 

In [16] we described in detail the procedure adopted to simulate the merger 
history of a cluster. We summarize here the basic points involved in this pro¬ 
cedure. 

In the PS formalism, the differential comoving number density of clusters with 
mass M at cosmic time t can be written as: 
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The rate at which clusters of mass M merge at a given time t is written as a 
function of t and of the hnal mass M' [20]: 
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where q is the present mean density of the universe, 5c{t) is the critical density 
contrast linearly extrapolated to the present time for a region that collapses 
at time f, and a{M) is the current rms density fluctuation smoothed over the 
mass scale M. For a{M) we use an approximate formula proposed in [21], 
normalized by assuming a bias parameter b = 0.9. We adopt the expression of 
6c{t) given in [22]. In this respect our approach is similar to that adopted in 
[23], 

In £g. 1 we plot a possible realization of the merger tree for a cluster with 
present mass of The history has been followed back in time up to red- 

shift z = 3. The big jumps in the cluster mass correspond to merger events, 
while smaller jumps correspond to what in the literature are known as ac¬ 
cretion events. During the merger of two clusters of galaxies, the baryonic 
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Fig. 1. Merger history of a cluster with present mass 10^® solar masses. The mass 
(y-axis) suffers major jumps in big merger events. Time is on the x-axis. 
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component, feeling the gravitational potential created mainly by the dark 
matter in the two cluster, is forced to move supersonically and shock waves 
are generated in the intracluster medium. 

In order to properly describe the physical properties of these shocks, we use 
here the approach introduced in [24,16]. We assume to have two dark matter 
halos, as completely virialized structures, at temperatures Ti and T 2 , and with 
masses Mi and M 2 (here the masses are the total masses, dominated by the 
dark matter components). The virial radius of each cluster can be written as 
follows 


^vir,i 


3M 


GMi 


^47r(1 T ^ 




( 3 ) 


where i = 1,2, pcr = 1-88 x cm“^ is the current value of the critical 

density of the universe, Zf^i is the redshift of formation of the halo i, Ac = 200 
is the density constrast for the formation of the halo and is the matter 
density fraction. In the right hand side of the equation we used the fact that 
Per = SHq/SttG, where Hq is the Hubble constant. The formation redshift 
Zf is on average a decreasing function of the mass, meaning that smaller 
structures form at larger redshifts, consistently with the hierarchical scenario 
of structure formation. There are intrinsic fluctuations in the value of 2 :/ at 
hxed mass, due to the stochastic nature of the merger tree. The formation 
redshift Zf is calculated here following [20] ^ . 


The relative velocity of the two merging structures, W, can be easily calculated 
from energy conservation: 


GM1M2 ^ ^ y‘2 _ GM1M2 

^vir,l T ^virg ‘2 2R\2 


( 4 ) 


where Mr = M 1 M 2 /(Mi + M 2 ) is the reduced mass and R 12 is the turnaround 
radius of the system, where the two subhalos are supposed to have zero relative 
velocity. In fact the hnal value of the relative velocity at the merger is quite 
insensitive to the exact initial condition of the two subclusters. In an Einstein- 
De Sitter cosmology this spatial scale equals twice the virial radius of the 
system. Therefore, using eq. (3), we get: 
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^ We adopt here as formation redshift the peak value of the distribution given in 

[ 20 ]. 
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where Zf is the formation redshift of the halo with mass Mi + M 2 . This ex¬ 
pression remains valid in approximate way also for other cosmological models 
[25]. The sound speed of the halo i is given by 

2 / 1 ^ 

Cs,i = MUa - ^)7r -, 

vir,i 

where we used the virial theorem to relate the gas temperature to the mass and 
virial radius. The adiabatic index of the gas is 7 ^ = 5/3. The Mach number 
of each cluster while moving in the potential of both clusters can be written 
as follows: 
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where 7 = M 2 /Mi < 1. The procedure illustrated above can be applied to 
a generic couple of merging halos, and in particular it can be applied to a 
generic merger event in the history of a cluster with hxed mass at the present 
time. 


The results of our calculations of 500 realizations of the merger history produce 
the Mach numbers plotted in £g. 2. The striking feature of this plot is that for 
major mergers, involving clusters with comparable masses (7 ~ 1), the Mach 
numbers of the shocks are of order unity. In other words the shocks are only 
moderately supersonic. In order to achieve Mach numbers of order of 3 — 4 
it is needed to consider mergers between clusters with very different masses 
iv 0.05). These events are the only ones that produce strong shocks, and this 
is of crucial importance for the acceleration of suprathermal particles, and, as 
discussed below, for the calculation of the spectrum of the diffuse gamma rays 
generated by the accelerated particles. 

Each merger here is assumed to be a two body event, namely the potential well 
is assumed to be dominated by the two merging structures. It may be argued 
that the merger between two objects may occur in the deeper gravitational po¬ 
tential created by nearby structures. In this case, the relative velocity between 
the two clusters, and also the related shock Mach numbers may be larger (or 
smaller) than those estimated above. Following [16] we briefly discuss here a 
simple argument that suggests that this problem should not be relevant for 
our purposes. We assume that the two clusters, with mass Mi and M 2 , are 
merging in a volume of average size Rgm where the smoothed overdensity is 
1 -|- (5 (5 = 0 corresponds to matter density equal to the mean value). Clearly 
the overdense region must contain more mass than that associated with the 
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Fig. 2. Distribution of the Mach numbers of merger related shocks as a function of 
the mass ratio of the merging subclusters. The upper strip is the distribution of Mach 
numbers in the smaller cluster, while the lower strip refers to the bigger cluster. 

two halos, therefore for a top-hat overdensity at = 0 we can write: 

+ S) = e(Mi + M 2 ), (7) 

where .^ > 1 is a measnre of the mass in the overdense region in excess of 
Ml -|- M 2 . In nnmbers, using Qm = 0.3, this condition becomes: 

(1 + 5) = 2eMl5i^^o^ (8) 


where M 15 is Mi -I-M 2 in units of 10^^ solar masses and Rgm = 10 Mpc i?io h 

If the clusters are affected by the potential well of an overdense region with 
total mass Mtot, the maximum relative speed that they can acquire is Vmax ~ 
2^GMtot/Rsm- Note that this would be the relative speed of the two clusters if 
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they merged at the center of the overdense region and with a head-on collision, 
therefore any other (more likely) conhguration would imply a relative velocity 
smaller than Vmax- In particular, the presence of the local overdensity might 
even cause the two merging clusters to slow down, rather than a larger relative 
velocity. In numbers 

Vmax = 1.1 X cm/s. 

Using the usual expression for the sound speed in a cluster with mass Mi we 
also get 

Cg = 8.8 X 10’^mI( 5 cm/s. 

Therefore the maximum Mach number that can be achieved in the i-th cluster 
is 




( 9 ) 


As stressed in the previous sections, the Mach numbers which may be relevant 
for particle acceleration are M. > 3, which implies the following condition on 

e > 5.8Mr5'M^(^i?io, (10) 

that, when introduced in eq. (8) gives: 

+ ( 11 ) 


Similar results may be obtained using the velocity distribution of dark matter 
halos as calculated in semi-analytical models [26] and transforming this dis¬ 
tribution into a pairwise velocity distribution, by adopting a suitable recipe. 

The probability to have an overdensity 1 -|- <5 in a region of size Rgm has the 
functional shape of a log-normal distribution, as calculated in [27]. Eq. (11) 
gives the overdensity l-|-(5 necessary for a cluster of mass M* to achieve a Mach 
number at least 3 in the collision with another cluster in the same overdense 
region. The probabilities P{5) as a function of the mass of the cluster Mj were 
estimated in [16]. For rich clusters, with masses larger than 5 x IO^^Mq (cor¬ 
responding to X-ray luminosities Lx > 4 x lO'^^ergf/s) the probability that the 
presence of a local overdensity may generate Mach numbers larger than 3 has 
been estimated to be less than 10“^, suggesting that our two body approxi¬ 
mation is reasonable, in particular for the massive clusters that are typically 
observed to have nonthermal activity. For smaller clusters, the probabilities 
become higher, indicating that the distribution of Mach numbers might have 
a larger spread compared with that illustrated in fig. 2. Note however that for 
small clusters, even the two body approximation gives relatively high Mach 



numbers, provided the merger occurs with a bigger cluster, simply as a result 
of a lower temperature and a correspondingly lower sound speed. 

Motivated by these arguments, in the following we keep the assumption of 
binary mergers. 


3 Shock acceleration during structure formation 


Merger related shocks may serve as cosmic ray accelerators. In this section we 
estimate the maximum energy achievable in one of these shocks for acceleration 
of electrons. We aim at using the most optimistic situation, to achieve the 
highest maximum energy for the accelerated particles. For this reason, we 
adopt a Bohm diffusion coefficient 

D{E) = -tlc = 3.3 X 10'^^E{GeV)B-^cmys, 

3 

where E is the particle energy in GeV and is the magnetic held at the 
shock in units of /j,G. Other possible choices for the diffusion coefficent are 
discussed in [28,16]. 


The acceleration time is dehned as follows: 


'^acc 



Di ^ D2 
-Ml U2- 


3D{E) r(r + 1) 
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where in the last step we neglected the jump in the magentic held at the 
shock, and we introduced the ratio r = ui/u 2 of the two velocities upstream 
and downstream the shock. The compression ratio r is related to the Mach 
number of the shock through the relation 

“ |Wf2 + 2‘ 

For relativistic electrons the main channel of energy losses is represented by 
ICS on the photons of the microwave background. The time scale for these 
losses is 

Tics ~ 4 x 10^®E(Ge'F)-^ s. 

Requiring that acceleration occurs faster than losses implies the following max¬ 
imum energy 

E^a. = 6.3 X l{)^B]l^h{r)-^/^us GeV, (13) 
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where h{r) = r(r + l)/(r — 1 ) and is the collision speed in nnits of 10 ® 
cm/s. It is worth stressing that larger diffnsion coefficients would signihcantly 
reduce the value of the maximum energy compared to that given in Eq. (13). 
This effect is discussed in [28] and [16]. 

It is instructive to evaluate the energy of the photons radiated by electrons 
with energy E^ax due to synchrotron emission and ICS respectively. The peak 
frequency for synchrotron emission is 

Umax = 1-4 X Blh{r)~^ul Hz. 

This corresponds to an energy 

hvmax ~ Wlh{r-)~^ul eC. 

The peak energy for ICS off the microwave photons is 

eics ~ TeV. 

The choice of a diffusion coefficient that is larger than the Bohm coefficient 
would considerably decrease the values of Umax and ejcs- The gamma ray 
production from ICS strongly relies upon the assumption of Bohm diffusion. 

The spectrum of the accelerated electrons is provided by the standard theory of 
hrst order shock acceleration, and in the linear regime is uniquely determined 
by the compression factor (or equivalently the Mach number) of the shock. In 
fact the spectrum can be written as a power law E~'^ with 7 = (r + 2 )/(r — 1 ). 
The slope tends asymptotically to 7 = 2 for strong shocks (r ~ 4). 


4 The diffuse gamma ray background from hierarchical clustering 


Electrons accelerated at each merger event lose energy through ICS on the 
photons of the cosmic microwave background, that are upscattered to higher 
energies, up to gamma ray energies. At each merger, two main shock sur¬ 
faces are generated [16], each one able to accelerate its own population of non 
thermal electrons. The balance between injection and energy losses drives the 
electrons toward their time independent equilibrium distribution, with a spec¬ 
trum one power steeper than the injection spectrum. The rate of gamma ray 
production from each merger, Q.y{E.y, Mi, M 2 ) (see [28] for details of the calcu¬ 
lation of from the electron spectrum) is thus the sum of the contributions 
from the two shocks, which in general depends upon the masses Mi and M 2 of 
the merging structures. Each shock generates a power law spectrum of gamma 
rays with two different slopes (determined by the compression factors at the 
shocks as explained in the previous section) and a cutoff determined by the 
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maximum energy of the accelerated electrons. The flux of gamma radiation 
(in units of phot cm“^s“^sr“^GeV“^) is then given by 

Zmax ^ OO 

0 ^ ^ 0 

M 

J dM' n{M, M + M', z)Q^{E^{l + z),M, M')AWr(M, M'), (14) 

0 


where S{z) = Y^r2m(l + zY + Ga, TZ{M, M+M', z) is the merger rate between 
clusters of masses M and M' at redshift 2 :, and Atmer{M, M') is the dnration 
of the merger between two clusters with given masses M and M', as dehned 
in [16]. 

In the literature the expression accretion events is often used [23,29] to describe 
the mergers of small subclusters with a large dark matter halo. This dehnition, 
implying an artificial separation between merger events and accretion, may be 
useful in other contexts, but here may be quite confusing. As discussed in §2, 
the so-called accretion events are the ones that generate high Mach numbers 
and therefore flat spectra of accelerated particles [16]. It would therefore be 
instructive but rather difficult to dehne the boundary between accretion and 
merger events. In the older literature, the concept of accretion onto a large 
scale structnre was discussed in detail (see for instance [30]) bnt it had a 
rather different meaning: a small perturbation in the density held grows as a 
result of gravitational instability, so that more matter falls onto the potential 
well. There is a radius, the so-called turnaround radius at which the inhow 
velocity is balanced by the expansion of the nniverse. An accretion shock is 
formed at a position that depends on time in the same way as the turnaround 
radius and as the virial radius of the structure. The shock surface carries the 
information of the virialization of the inner region of the cluster. This type 
of accretion is also called secondary infall, meaning that matter accretes on 
a potential well which has already been formed. In the following, we adopt 
this simple approach and calculate the gamma ray production due to particle 
acceleration at the accretion shock, and compare it with the result of the 
gamma ray production from mergers. 

For simplicity let us assume that the shock is located exactly at the virial 
radius Rmr- The total energy per nnit time flowing across the shock is then: 

1 

Ltot = -Pcr-^(1 + zfvjf47rRli^, 

Z i Lyyi 

where Pcr^m is the matter density of the nniverse, is the baryon fraction 
and Vff is the free-fall velocity of the matter at the distance of the shock radius. 


11 



The secondary infall just described is a simplification of the mass flow through 
large scale shocks in the hlamentary structures seen in N-body simulations. 
Although the geometry is different, the total energy crossing the shock per 
unit time and per unit surface should not be very different from the same 
quantity calculated for spherical inflow. 

The accretion shock, by dehnition, propagates in a cold (non-virialized) medium, 
therefore its Mach number may be very high, although the typical speed of 
the shock is of the same order of magnitude as the merger speed of two clus¬ 
ters. The exact value of the Mach number depends on the temperature of the 
medium before entering the overdense virialized region. If we take for such 
temperature the range of values T = 10^ — 10® K, and a typical shock speed 
of ~ 10® cm/s, the corresponding Mach numbers range between 10 and 100. 
These Mach numbers, being much larger than unity, correspond to spectra of 
accelerated particles which are ~ E~‘^, with a cutoff at the maximum energy 
(Eq. 13). If the intergalactic medium were pre-heated before the gravitational 
collapse, these Mach numbers could be lower [31]. 

The rate of gamma ray production from the accelerated electrons, M, z), 

can be calculated in the usual way [28]. Given Q^{E^, M, z) (in phot s“^ GeV“^) 
in an accreting cluster of mass M at redshift ; 2 , the diffuse flux is 

Zmax ^ OO 

° 0 ^ ^ 0 


The diffuse flux of gamma radiation from mergers (dashed line) and from 
accretion (dotted line) is plotted in Fig. 3, where an acceleration efficiency 
(for electrons) of 5% is assumed. The observed EDGRB is the shaded region 
[1]. In the same hgure we plot for comparison the predictions of Ref. [14] (solid 
line), where our same acceleration efficiency was adopted. The meaning of the 
dash-dotted line will be explained below. 

Three conclusions are evident: 

1) the flux of gamma radiation from both mergers and accretion is a factor 
~ 10 smaller than the observed EDGRB and smaller than the flux predicted 
in [14,15], by the same factor. An acceleration efficiency of the order of 50% 
should be adopted in order to reproduce observations. This would be unrea¬ 
sonable for electrons as accelerated particles, and would violate our initial 
assumption of shock acceleration in the linear regime (no backreaction of the 
accelerated particles on the shock); 

2) the gamma ray diffuse flux from mergers is at the same level as that due 
to accretion (secondary infall); 
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Fig. 3. Diffuse gamma ray emission from structure formation. The shaded area is 
the result EGRET observations. The dashed line is the result of our calculations 
for mergers while the dotted line is the flux of gamma rays from accretion. The 
dash-dotted line assumes a minimum mass of the merging halos o/lO^^M©. 

3) all predicted spectra are approximately power laws with slopes between 2 
and 2.1. 

The flatness of the spectrnm may give ns the key to nnderstand the differences 
between onr resnlt and that of [14,15]. The anthors of Ref. [14,15] assnme that 
all the shocks are strong, so that the spectrum of the accelerated particles is 
fixed to E~'^. In our approach the Mach numbers of the shocks, as well as 
the energy flux through each shock are calculated self-consistently, so that the 
spectrum of the diffuse radiation is the result of the superposition of spectra 
with different slopes. However, although major mergers are energetically very 
powerful, they generate steep spectra of accelerated particles [16]. Therefore 
the contribution of relativistic electrons from these mergers is small compared 
with that of smaller, less energetic mergers, which however produce flatter 
spectra. This is the reason why the gamma ray spectra resulting from mergers 
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have almost the same spectrum as that predicted in [14,15]. Our absolute 
normalization is however a factor ~ 10 lower, which may suggest that merger 
related shocks are not always strong, and in fact they are almost always weak, 
as illustrated in Fig. 2. 

Our results can be better clarihed by using Fig. 4. The upper panel shows 
the average normalized energy flux per unit time through the merger related 
shocks of a cluster with mass IO^^Mq (solid line), IO^^Mq (dotted line) and 
IO^^Mq (dashed line) at redshift z = 0, as a function of the mass of the 
merging subcluster (here M'/M < 1 is the ratio between the masses of the 
two subclusters). The curves represent the energy flux contributed by mergers 
with mass ratio larger than M'/M. The energy flux sums up to ~ 70 — 80% 
of the total for subcluster masses larger than ~ O.IM. This implies that the 
energy flux that crosses the shocks formed during mergers of the cluster with 
mass M and subclusters with masses smaller than O.IM is small (20 — 30%). 
In other words, the energy flux is dominated by major mergers. 

In the second panel of Fig. 4 the energy flux is plotted for a cluster of hxed 
mass of IO^^Mq at three redshifts, z = 0 (solid line), 2 : = 0.5 (dotted line) and 
z = 1 (dashed line). The same conclusions explained above hold here. 

The third panel is the most interesting: it represents the normalized energy 
flnx throngh the merger related shocks that contribnte to the diffuse gamma 
ray background above 100 MeV, for clusters of masses as labelled in the npper 
panel, at 2 : = 0. It is immediately clear that most of the contribntion to the 
gamma ray emission is provided by mergers with small mass ratios, M'/M < 
10“^, namely the ones having the larger Mach numbers (see Fig. 2). 

Summarizing, while most of the energy flows through shocks associated to 
major mergers, the energy flux that contributes to the gamma ray background 
is the one that crosses strong shocks, occurring when a large cluster encounters 
a snbcluster with ~ 0.01 times the mass of the larger clnster. This may explain 
why the diffuse gamma ray background as derived in the present paper is 
substantially lower than that estimated in previous calculations [14,15]. In a 
recent paper [32], a reevalnation of the diffnse gamma ray backgronnd from 
large scale strnctnre shocks was carried ont and there seems to be there a 
closer agreement with the conclusions of onr calcnlations. In [32] many issnes 
were discussed as possible reasons for the discrepancy with the results of [14]: 
one of the points that the anthors correctly hnd ont is that simnlations do 
suggest the formation of weak shocks, althongh difficult to identify. Another 
nnmerical calculation was also carried ont in Ref. [33]. In [32] the anthors 
emphasize the difficulties in the identihcation of shocks with Mach nnmber 
below 10 (and the impossibility to detect shocks with Mach nnmbers below 
~ 3 — 4). The Mach nnmber distribntion obtained in [32] presents an artihcial 
peak at the threshold of detectability (A4 ~ 4). A peak is seen also in [34], 
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but it is argued that it is not an artifact of the numerical procedure. In our 
semi-analytical approach no peak is found at ~ 4, while the Mach number 
distribution seems peaked at ~ 1.5 116|. 

Another ingredient introduced in our calculation (as well as in [32]) but not 
in [14] is the redshift dependence of the 7 —ray emissivity. This also induces 
7 —ray diffuse fluxes smaller than those in [14], as also pointed out in [32]. 
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Fig. 4. Upper panel: normalized energy flux per unit time through the merger related 
shocks of a cluster with mass 10 ^^Mq (solid line), 10 ^^Mq (dotted line) and IO^^Mq 
(dashed line) at redshift z = 0 with clusters with mass larger than M'/M. Middle 
panel: Same as above for a cluster of mass at redshifts z = 0 (solid line), 

z = 0.5 (dotted line) and z = 1 (dashed line). Lower panel: normalized energy 
flux through the merger related shocks that contribute to the diffuse gamma ray 
background above 100 MeV, for the same halos as in the upper panel. 

The flux of diffuse gamma rays due to accretion (secondary infall) has a spec¬ 
trum which is exactly E~‘^ because the Mach number of the accretion shock 
is always much larger than unity. It is somewhat surprising that the diffuse 
gamma ray background contributed by electrons accelerated at the accre- 
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tion shock is comparable with that produced in merger events (the latter, 
as stressed above, is dominated by mergers between clusters with M'/M <C 1, 
that in some literature are indeed dehned as accretion events). In this respect 
however some additional discussion is required: for a cluster with mass 
a mass ratio M'/M ~ 10“^ corresponds to a substructure with mass 
comparable with the mass of our galaxy. This clearly does not make physical 
sense. Galaxies move within the intracluster medium without their medium 
being shocked. It is more likely that a bow shock is formed in front of the 
galaxy, due to the internal pressure of the galactic medium [35]. Actually sim¬ 
ulations show that large galaxies penetrating the intracluster medium of a rich 
cluster can even be completely stripped of their gas content [35]. This suggests 
that a low mass cutoff should be imposed in the calculation of the gamma ray 
diffuse background from cluster mergers. The dash-dotted line in Fig. 3 has 
been obtained by considering only structures with virial masses larger than 
corresponding to galaxy groups. The diffuse background of gamma 
rays is a factor ~ 10 smaller in this case and is slightly steeper in spectrum. 
This happens because the main contribution comes from mergers between clus¬ 
ters with masses and M = ~ Clusters 

with masses as large as are already on the tail of the Press-Shechter 

distribution even at z = 0, therefore the corresponding contribution to the 
diffuse background is suppressed. On this basis, the dash-dotted line in Fig. 
3 is the most realistic estimate of merger shocks to the diffnse gamma ray 
background, amounting to ~ 1% of the observed EDGRB (this result agrees 
with the estimate in [36]). On the other hand the strong shocks associated to 
accretion of matter onto a cluster may generate a gamma ray background as 
large as that plotted as a dotted line in Fig. 3, and this contribntion remains 
at the level of ~ 10% of the observed EDGRB. This may be considered as 
the most realistic prediction of the contribution of clusters of galaxies to the 
EDGRB. 

Our results are in agreement with the recent estimate of the diffuse gamma ray 
flux from rich clusters carried out by cross correlating high galactic latitude 
EGRET data with the location of Abell clusters [37]. 


5 Conclusions 


We calculated the contribution of structnre formation to the EDGRB radia¬ 
tion. We hnd that this contribution is only 1 — 10% of the observed flux of 
the alleged extragalactic radiation above 100 MeV as measured by EGRET. 
The calculation has been carried ont in two different scenarios, one that relies 
upon the hierarchical scenario for structure formation, and the other that is 
based on the secondary infall (accretion) of matter onto a potential well which 
has already been formed. 
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In the hierarchical approach, large structures are formed by merging of smaller 
halos whose mass is dominated by dark matter. The baryon components of 
these halos, moving supersonically, develop shock surfaces that in principle can 
accelerate particles to relativistic energies. We evaluate the Mach numbers of 
these shocks following the recipe introduced in [16,24], that allows us to self- 
consistently calculate the spectra of accelerated particles at each merger event. 
This distinguishes our approach from previous calculations [14,15] where the 
shocks were all assumed strong (inhnite Mach number), so that the spectra 
of accelerated particles were by dehnition E~‘^. We hnd that this assumption 
may lead to incorrect conclusions, as evidenced by the comparison between 
our results and those in [15]. 

In [16] we investigated the role of protons first accelerated and then diffusively 
conhned in large scale structures [12]. We find that the spectral shape of the 
protons is very steep and is unlikely to produce a relevant effect on high energy 
radiation generated by clusters of galaxies. This conclusion is mostly due to 
the fact that major mergers, that energetically dominate over smaller merger 
events, generate weak shocks and therefore steep particle spectra. 

In the present paper we focused our attention on electrons as accelerated 
particles. Their ICS energy losses were in fact proposed [14,15] as responsible 
for upscattering the photons of the cosmic microwave background to gamma 
ray energies, therefore generating a diffuse background of gamma radiation 
accompanying the process of structure formation. While previous calculations 
suggest that the observed EDGRB may be saturated by the contribution of 
particles accelerated during structure formation, we hnd here that at most 
10% of the observed background can in fact be explained in this way. More 
recent numerical calculations [32,33] seem to lower previous predictions. 

Although structure formation is generally believed to follow the hierarchical 
picture outlined above, some secondary infall of matter onto forming struc¬ 
tures must occur, and is in fact observed in numerical N-body simulations in 
the form of hlamentary-like accretion hows. In order to account for this con¬ 
tribution we adopt a simple model, similar to that proposed in [30], in which 
an accretion shock is formed at approximately the virial radius of a struc¬ 
ture that is accumulating matter from the expanding universe. This shock, 
as those formed in the hlaments, propagates in the cold unshocked medium 
and may have very large Mach numbers, of order 10-100 (possibly lower if 
pre-heating occurs [31]). Particles accelerated at this type of shocks have the 
hattest spectrum allowed in the linear theory of shock acceleration, namely 
E~‘^. The dihuse gamma ray background due to ICS of electrons accelerated 
at accretion shocks is also of order 10% of the observed EDGRB, and may 
be larger that that due to mergers. The reason for this result is the following: 
although the energy flux through merger shocks is larger than that crossing 
accretion shocks, the former is mainly converted into very steep spectra (due 
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to the weakness of the shocks), while the latter is more easily channelled into 
particles that may contribute to gamma rays with energy above ~ 100 MeV, 
because of the flatter spectra. 

In the perspective of future work in this held, it seems to us that priority 
should be given to improve N-body simulations in order to have a better 
handling of the shocks with intermediate strength generated during structure 
formation. This will allow a self-consistent treatment of both the gas heating 
and the acceleration of particles at these shocks, and make a solid case in favor 
or against clusters of galaxies as sources of high energy gamma radiation. 

The results presented in this paper suggest that clusters of galaxies, and more 
in general structure formation, do not contribute appreciably to the EDGRB. 
However, several clusters could be observed as single gamma ray sources by 
future experiments such as GLAST (Blasi and Gabici, in preparation) and 
provide useful information on the non thermal history of these large scale 
structures. 

Acknowledgments: We are grateful to U. Keshet, G. Dermer and F. Stecker 
for useful comments. 
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